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As a possible approach to the neutrino oscillation on the basis of quantum field theory, the 
expectation values of the flavor-neutrino currents are investigated by employing the finite-time 
transition matrix in the interaction representation. Such expectation values give us in the simplest 
form a possible way of treating the neutrino oscillation without recourse to any one flavor-neutrino 
states. The present paper is devoted to presenting the formulation and the main structures of the 
relevant expectation values. 

■ PACS numbers: 14.60.Pq 

o : 
o ■ 

■ I. INTRODUCTION 
^\ 

D ' Many papers have been published concerning the theoretical basis of the standard formulas of the neutrino oscil- 
lations. As to the works appeared till the middle of the year 2001, we can find the list of references in the review 
\Q ' article written by Buete[3|; in this article he gave his view by arranging the controvercial points in the neutrino- and 
^SJ . meson-oscillations. 

T— I ■ 

^ . A. Proposal of problems 

I> ■ 

■ The analyses of the neutrino experimental data have been usually done on the basis of the standard formula of 
CnI \ the neutrino oscillation. In the present work we concentrate on the following two problems included in the ordinary 

■ derivation of the standard formula of neutrino oscillations in vacuum 0, Q . 

"nI ' First we summarize the ordinary derivation of the standard formula in vacuum0,l3. Let's consider the neutrino with 
momentum k and helicity r = — 1 , produced through the charged-current weak interaction. Such a flavor-neutrino 
state is assumedto be a superposition of the mass eigenstates, expressed as 

Oh. 

(i..(k)| -^(;.,(k)|i.,(k))(i.,(k)|, a = (e,/i,...),j = (l,2,---); (1.1) 



Oh! 

L| ' the matrix U — [Uaj] with U^j = (i'o.(k)|i'j (k)) is unitary. 
, By employing the time evolution of (i/j(k)| , 

^, (^,(k;0| =e--^('^)*(^,(k)|, ..,(k) = y^k2 + m2, (1.2) 

the flavor-neutrino state is seen to evolve in time as 

(^p(k;t)| =^C/„e— ^('')*f/;,.(z..(k)| 

= Y,Axp{\^;t){vx{\^)\, (1.3) 

A 

where 

Xp(k;t) = (t/e-'''^<^-«('')*C/t) , 
/^i(k) 

^diag = 

V 



c^2(k) 



(1.4) 
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Then, the transition probabihty is 



P,^p(i) = |Ap(k;Ol' 



J2 Uaj U*„ [/^j (i^) C'))* . 



Under the condition of the neutrino velocities 



1 (the hght velocity in vacuum), 



(1.5) 



(1.6) 



the time t which is regarded as the time of neutrino propagation is replaced by the distance L from source to detector 
approximately. Then we have 



(1.7) 



here we used LOj (k) — uji (k) 



2|k| 



2|k| 



, which holds for extremely relativistic neutrinos. 



M 



Now we mention the two problems included in the derivation described above. 
Problem I: The matrix U = [C/o-p] is defined as the matrix which diagonalizes the mass matrix in the flavor bases 



M = 



ee "^e^ 



fie \ , M ^ M'^ 



V 



' mi 



M, 



diag 



(1.8) 



(1.9) 



leads to 



t/M|k| + 



|k| 



'"diag 

2|k| 



^diag(k). 



(1.10) 



in the extremely relativistic case. 

The important point is that the matrix U is set equal to the mixing matrix Z^^'^ = \z'J^\ between the two kinds of 
neutrino fields, i.e. 



for 



vf{x) 



( Veixy 



Vi(.t)^ 



(1.11) 



(1.12) 



where Z^/^ is defined by 



diagi 



M is the mass matrix for the "free part" Co{x) = —vf{x){^ + M)vp{x) in the total Lagrangin 

C(x) = Coix) +Cint{x). 



(1.13) 



(1.14) 



Here, £i„t(x) is assumed to have no bilinear term and no derivative of the neutrino fields. The Hamiltonian Ho{x) 
obtained from Cq{x) is 



no{x) = Df{x){'j ■ d + M)vf{x) 

^Yvj{x){~^ ■ d +mj)vj{x), 



(1.15) 
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which means the energy is diagonal with respect to ^'j(a;)'s, irrespectively of the momentum. Thus, in spite of the 
situation that , apart from a trivial phase factor, we can take 

there is a certain inconsistency. It is often asserted j4| that such an inconsistency disappears in extremely relativistic 
neutrinos as in existing experiments ; therefore, the state with U = Z^/"^ reflects correctly the present exper- 
imental situation. It is , however, not so simple conceptually for us to solve such an inconsistency from the field 
theoretical viewpoint as noted in Refs.|^-0. The problem with which we are now confronted is what relation the 
state has with the field-operator relation; in other words, what definition should be given to the one-particle state of 
flavor neutrino properly under the condition vf{x) = Z^^'^vm{x). 

In connection with this problem, various authors have proposed [9l|-|lll| that the neutrino oscillation shoud be 
investigated on the basis of the quantum field theory by examining transition amplitudes, in which flavor neutrinos 
appear only intermediate states between source- and detection-interactions, as schematically shown in FIG^ In such 
approaches, it is unnecessary for us to prepare any one-particle state of flavor neutrino from the beginning ; therefore, 
on the contrary to the standard formulation, we are free from the trouble concerning Problem I. 




FIG. 1: Real experimental situation 



Problem II : The replacement of ct by the source-detector distance L for obtaining (1.7) reflects such a situation 
that Pa-^p{L) mimics the usual experimental circumstances more closely than Pa-^p{t), because the times of reactions 
are not known and we measure distances rather than times. This simple replacement causes some doubts concerning 
the consistence of this procedure with the approximation Wj(k) — Wi(k) ~ ATO|j/(2|k|) . 

Concerning a proper definition of space-dependent probability Pa^p{L), an interesting approach called "the current- 
density approach" has been proposed by Ancochea et al.|l2l| ; some modification is found in Ref.|l3|. In this quantum 
mechanical approach, the oscillation formulas P^^p(L), as a function of the source-detector distance L, are obtained 
after integrating the probability current density over a surface dA arround a detector as well as the time duration of 
measurements ; 

Pa^p{L)= f'dt f dS.j(x,t). (1.17) 

Jti JdA 

By employing the appropriate probability current densities of K'^ and K'^ constructed with the use of wave-packet 
functions of and K^, the usual formulas have been shown to be derived. 

Blasone et al. 0| considered the neutrino oscillation along the line of the current density approach. The current 
density jp(x, t) is replaced by the expectation values 

(z.,(x/,t7)|jp(x,i) \v^{^i,ti)), (1.18) 

where |z^ct(x/, i/)) is an initial wave-packet state of the flavor constructed on the flavor vacuum in accordance with 
the formalism given by Blasone and Vitiello 0, ■ It is pointed out that in low-energy region of neutrinos certain 
significant deviations from the standard formula may appear. 



B. Purpose and starting formula 

The purpose of the present work is to investigate the neutrino oscillation by taking into consideration the two 
kinds of approaches explained in connection with Problems I and II. We consider the formula (|1.17|) with jp(x, a;°) 
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replaced by a certain kind of expectation value in the quantum field theory. In order to do this, we notice that in the 
interaction representation the expectation value of a physical observable F{x) at a sp ace-time point x = (x, x*') with 
respect to a state is given, in accordance with e.g. Umezawa's textbook by 



Fix) = {'f{x"j)\S-\x",x"j)F{x)S{x°,x'j) 



(1.19) 



where 



\^ix'))=S{x",x',)\^ix",)), 



m=0 



A2 



d'^yrnHint{yi)Hint{y2) ' ' ' Hint{y„ 



(1.20) 



We consider the simplest case, where one 7r+ decays into z^^ and /I; for convenience we write tt"*" v„ -\- l^. As 
the initial state [^^(x")), we adopt a wave-packet state which is to be explained later in Sec. III. So we examine 
the expectation value of the flavor-neutrino current j^^j{x), which is to be defined below, with respect to the state 




(L»(v„)/(r.>) 



FIG. 2: 



The content of the first-step investigation is schematically shown in FIG [3 Thus this is seen to give the simplest 
case, in which the oscillation probabilities are calculated in accordance with the quantum field theory as well as 
without recourse to any flavor neutrino state. After a small modification, it will be seen that we can take into account 
the observation of ^+ as schematically shown in FIG|31 




FIG. 3: 



The relevant weak interaction in (|1.14|l is written as 

Cint{x) = ^ [vp{x)Jp{x) + Jp{x)vp{a 



(1.21) 



Jp{x) = 5/5v{x) ■ Ant = if-„pV°-lp{x)da(j){x), 
Jp{x) = dCint/Sv{x) = if* lp{x)v°-da(j)\x), 



(1.22) 
(1.23) 
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with v"' = 7°(l-l-75)- (We use the Kramers representation of 7-niatrices as explained in Refs.|^ and 0) Following 
the natural definition of " currents" [l7j , we consider the infinitesimal phase transformation 

M^) ^ ^^f(^) - e^'^^'^Kpix), (1.24) 

with 

(ipe{x) 

and obtain 

£(2;) ^ C\x) ^ C{x)+6£{x), 
6C{x) = ^wf{x) ii^davix) + [M, ^{x)]) vf{x) + 5Ci,,f (1.25) 

By remembering Ci-at has no derivative of J'i?(x)-field, we obtain the flavor neutrino current as 

dSCjx) 
d{daipp{x)) 
-iDp{x)-f°'i'p{x) 

-^Y,zlfzli;'D,{x)r^l[x)■ (1.26) 
0,1 



its 4-divergence is given by 



^ -a / N d5C{x) 
^"^('')^^^-5(^ 

= -z [{DF{x)M)pVp{x) - Pp(x)(Mj.i.(x))p] + ^^1^. (1.27) 

The right-hand is called the " diffusion term " [T^ Il3|. In the interaction representation, in which the j^M(a;)-field 
satisfies the free equation 

(^ + Afdiag);^A./(x) =0, (1.28) 
the total flavor current of neutrinos is seen to be conserved; 



p 



dajl)=0. (1.29) 



The remaining part of the present paper is organized as follows. In Sec. II. we perform some model calculations to 
see the relation of the current-operator expectation values (1-19) to corresponding transition amplitudes. In Sec. Ill, 
we investigate the structures of the expectation values for j^^-^ (x) with respect to the wave-packet state, decaying 

into and /o-- In Sec. IV, we discuss the structures of the same kind of the expectation values in the boson model 
for simplicity, in which the observation of a charged boson in decay products is taken into consideration, and in 
Sec.V summarize results and future tasks. The detailed analyses including numerical calculations are to be left to the 
subsequent paper 0|. 

II. RELATION OF CURRENT EXPECTATION VALUE TO TRANSITION MATRIX ELEMENT 

A. The case with no mixing 



We perform the simplest model calculation with the aim of seeing the relation of a current expectation value to the 
absolute square of the relevant transition matrix element. First we examine the case of no mixing interaction. 
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Let's consider the simple situation where (pseud-) scalar particles A, B, C have the weak interaction 

Cint = fABC(l>c{x)(f>B{x)(j)A{x) + H.C. 

= -Hi„t(a;). (2.1) 

An yl-particle described by the complex field (j)A{x) decays into B and C particles. We write the current of the C-field 
as 

= i^l,^i,^c{x). (2.2) 

We execute, in the interaction representation, perturbative calculations in the lowest order of Hint by employing 
the plane wave expansion 



1 



, f^A,B,C; (2.3) 



here, px — p -x — Ef{p)x^ , ^/(p) — i/p^ + m^; the commutation relations among the expansion-coefficient operators 
are 

[a/(p),at(q)] = [/3^(p), /^t (q)] ^ 2Ef{p)Sfg6{p - q), others = 0; (2.4) 
the vacuum |0) is defined by 

a<,(p)|0) =/33(p)|0) =0 for 9^A,B,C. (2.5) 
the number operator of ^p-field is given by 

Ncix°)^i: fd'xjPix) : 



d'k^^^^[al,ik)acik) - plik)Pcm. (2.6) 

(The symbol : : means the normal ordering.) 

In the following we examine the forms of expectation values of the current and the number operator with respect 
to one j4-particle state 

|^(p;x°)) = 5(x°,x?)at^(p)|0). (2.7) 
We define the expectation value of the current with respect to the state (|2.7ll as 

f(A(a;?;C-current(a;))b = (0| a^(p)5(x", x?)"! : j(=)(a;) : 5(xO, x?)at^(p) |0)™-, (2.8) 

and also that of the number operator as 

{Nc{x^))^^p^^o^ = i Jd^x£{A{x"j;C-cuTrentix))i. (2.9) 

We take out only the connected part in the expectation values, as designated in R.H.S. of 12.8|l . In the lowest order 
( i.e. the second order ) of the weak interaction, we have 



(iVc(x°)) =(0|«^(p) / d^zHUz)Nc{x'') Ailint(2;)ali(p)|0), (2.10) 
which consists of two parts ; 

{^C{X')) Ai,,.»,) = ("C.(-°))^(,,,0) + (-C(a^°))^(,,,0) , (2.11) 
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where 

nc 



Jd'k^^L_aU^)ac{k), nc^- Jd'k^^^^plik)(3c(k). (2.12) 



The first part, represented as FIG^(a), is given by 

(^^(^°)>A(p,.o) = £ d^^£ d^y\I^^e-^P^+wy(^Q\ ^U^)My) |0>(0| M^MUv) |0> 



f f dz°dy°V\fABC? fd'k 
Jo Jo J 



»(£;^(p)-£;B(q)-£;c(k))(^»-a») 



(2^)62£;s(q)2^c(k) 



(2.13) 

q=p k 



here, t = — x'j, V = Jd^xl. Similarly, the second part in R.H.S. of H2.11|l . represented as FIG0](b), is given by 

i^cix'')) 0.= d'z d'y\fABC?{Q\Mz)<t>Uv)n 

X (0| <j,l{z)nc^c{y) |0)(0| a^(p)0^(z)</>t^(y)at^(p) |0)™" 

= - f fdz'dy^V\fABC? fd'k' 
Jo JO J 



rt ft r p-i{EA{p)+EB{q)+Ec{-k.mz°-v°) 
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,0| ,|.,0 



(27r)62£;s(q)2£;c(k) 



(2.14) 



q= p k 

By adding the damping factor e"^'^"'^''^' ' , e > 0, to the integrands of 1)2. 13|) and (|2.14l) . we obtain 



{2T:f2EB{f\)2Ec{\^) 

gi(£;A(p)--Ef3(q)--Ec(k)+ie)t _ ^ g-i(BA (p)-Bf3 (q)--Ec (k) -ie)i _ ^ 



X 



Ea (p) - Eb (q) - Ec (k) + ie Ea (p) - Eb (q) - Ec (k) - le 



(2.15) 



q=p-k 



("C(X°))^, 0, = -V\fABC? Id'k ^ 



A{p,xO) 'u/i^oiy- (27r)62£;B(q)2£;c(k) 

^-i{EA{p)+EB{<i)+Ec{]i)-ie)t _ I gi(SA(p)+BB(q)+Bc(k)+ic)t _ ^ 



X 



Ea{p) + Eb (q) + Ec (k) - (p) + (q) + Ec (k) + ie 



(2.16) 



q=-p-k 

Thus we see that due to the higher configuration in the intermediate state of FIG^l (b) in comparison with FIG01 



(a), the main contribution to {Nc{x°))^. a^ comes from {nc{x°)) ^, ,^0)', 



I.e. 



'A(p,x5) ^-""^ x'-ov-- W A{p,x°)' 

(^c(x"))^(^,,o)-(nc(-°))^(,,,o). (2.17) 

On the other hand, the transition amplitude from the initial A state at x^ to B{q) and C{k) at x^ is , in the first 
order of the weak interaction, given by 

M{A{p,x°)-.B{q)+C;x°) 

r" 

= (0| PB{<l)ac{k)i / d^zfABC^U^)^B{z)Mz)a\{p) |0) 



/•x° iz{p~q~k) 
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from which we obtain 



J(p-q-k) 



'(2^)62£;s(q)2^c(k) 



'(2;j0^yOgi(z°-y«)(£;A(p)-BB(q)--Bc(k)) 

Diagramatically we may express the above equaUty as FIG[S1 



(2.19) 



,0 „„0 



^11 




2EB(ci) 2Bc(k) 



FIG. 5; Relation of the expectation value of number operator to the transition probability. 

A derivation which is somewhat more general than that explained above is given when we assume the fic- 
contribution to the expectation values p. 8(1 is smaller than nc-contribution, because the former arises from higher 
configurations ; i.e. 

^0 ^On t 



(iVc(a;°))^(p,,o) = (0|a^(p)5-i(x",a;?)[nc + nc]^(x°,xO)a^(p)|0)^°", 



2Ecik) 

d^k 
2-Bc(k) 



(0| aA(p)5-i (k) |0)™" (0| ac (k)5(a:° , (p) |0)^ 



^ |(C(fc);X|5(x°,x?)|aA(p)) 



con 1 2 



(2.20) 



phase vol 



A similar consideration applies also to the weak interaction process, in which a source and a target are included 
and the particle C scattered from the target is observed, as shown in FIGlHl: 



E 

phase vol 




phase vol 



X 



A 



,0^ C B 



FIG. 6: Relation of the expectation value of number operator to the transition probability. Both include production and 
detection processes. 

In the case where X is one D-particle state, RHS of FIG El corresponds to the process examined by Yabuki and 
Ishikawa jT]| when the mixing exists in the intermediate state. 
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B. The case with mixing 



Along the aim of the present paper, it is important for us to consider the case where there exists a "flavor" degree 
of freedom with mixing among (pseudo-)scalar fields <j)cp{^)i P = 1; 2, ■ • ■ . The relevant Lagrangian is 



£(x) = Co{(l)A- and 0s-fields) - (SVcl^^) ' 5ft0c(a;) - (p]j{x)M^<j)c{x)) 
+ X! {fABCp(t)cp{x)(jjBix)(l)Aix) + iJ.c.j , 







/ 


c{x) = 


4>p{x) 






\ ■■ 1 


V 



{My = 



(2.21) 
(2.22) 



Similarly to lfrTT|) - ((T~T^ . we have 



0c(x) = zV20<,(a:)W, 4^'^)(x) 



(t>C2{x) 



7l/2^ rrl/2 



(2.23) 



(2.24) 



We assume m| > for all j = 1, 2, • • • , and now consider the interaction representation, where 0g-fields, g = A, B, C^-'s 
, satisfy the free equations of motion with masses m^'s, and Fock space is constructed on the vacuum |0), defined in 
the same way as (|2.5|) . 

We examine the structures of the expectation value of Ncp{x^) defined by 



Ncp{x^)^i: d?xjf''\x) : 



with 



P = e, //,••• ; 



(2.25) 



(2.26) 



here, needless to explain, the symbol : : represents the normal ordering with respect to 
{acj(k), /3c'j(k), j = 1, 2, • • • } and their Hermitean conjugates. Concretely we have 



Ncp{x-)=Y.zTKNd'k\ 

3,1 



1 1 

+ 



X [«t^^.(k)acKk)e^(^''^('')-^^'('^»^" - /3t,^(k)/3c,(k)e^(^^'('')-^«^-(''))-"' 



+ 



i— ) /3c, {-k)aci (k)e-(^-^ ('^)+^- C^))-" 



(2.27) 



here, /?cj(q) with {q = — k, qo = ^^^(k)} is written simply as (3cj{—^)- 
Similarly to (|2.8|) , we write the expectation value of jl'^^'' (x) as 



£(A(x?);Cp-cur(a;))fc ^ (0|aA(p)5(a;°,x?)-i : '''(x) : 5(a;0,x?)at^(p) |0)™^ 



(2.28) 



which leads to 



(iVcp(a;°))^(p,,o) ^ * / d'x£{A{xy,Cp-cmix)U. 



(2.29) 
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As we see from (|2.27ll . Ncp{x'^) consists of 4 terms, which we write Nc p{x'^)^"\ n ~ 1,2, 3, 4, in order; 

1 1 



„0n(3) 



,0n(4) 



(2.30) 
(2.31) 



In the lowest (i.e. the second ) order contribution with respect to the weak interaction in (|2.21|) . we have for n — 1,2, 



(2.32) 



(7Vcp(x")^ Oa(p,.9) = (0|«a(p) / d''zHUz)Ncp{xT' / d%HUyMi>)W 

' Jx° Jx° 

(iVcp(x")^'+'%(p.,o)=z2(0|a^(p) / / d''zdS[H,^t{z)H,MNcp{x"r 
+ Ncp{x"f ^H,MHintiz)]a\{p) |0)™". 



^^int(a;) = -fABCa<l>Ccr(^)'t'Bix)4>A{x) +H.C.. 

Then we obtain, as the respective contributions from the diagrams FIG's|Z|(l) and (2), 



X e 



3.1 

i(£;cj(k)-BCT(k))t 



dz^'dy'' 



^0 



(2^ 



2EBiq) 



X gi(BA(p)-BB (q)--Ec3 ik))z° ^^(-Ea (p)+Eb (q)+Bci (k))a° 



{Ncp{x''r)A(,..x^l)^-V\fABC.?Y.Ki Ki^ifz^ 

hi 



q=p-k 

d^k 



^0 



{2n) 



4 



i?c,(k) i?ci(k) 



2EB{ci) 



X g-«(BA(p)+-EB(q)+-Ec3(k))z''gj(iJA(p)+BB(q)+Bci(k))y« 



q=-p-k 



here, t — 



(2.33) 



In the following, we assume only one flavor cr is dominant in Hintix), similarly to n'^ decay (i.e. 7r+ /i+ + ^^/j) j so 
that we take 



(2.34) 



(2.35) 



(2.36) 



(1) 




(2) 




FIG. 7: Diagram representing II2.35|I and l|2.36^ . respectively 



Due to the same reason as explained in the preceding subsection A, (A^cp(2;*')^^'')/i(p,x5) contributes dominantly 
over (iVcp(a;°)^^^)A(p,:EO)- 
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As to {Ncp{x'^) )a{p.x°)^ the daigrams represented by FIG |H1 contribute to it, the form of which is concretely given 

, , Jo JO 



dpk 1 

(27r)6 4 



i?c,(k) Eci{\^) 



1 



„i(£;cj(k)+Sc!(k))t 



2i;B(q) 



gj(-BA(p)-BB(q)-Bci(k))2:''gi(£;A(p)+-BB(q)--Ecj(k))y" 



q=p-k 



q=-p+k 



(2.37) 



(1) 




(2) 




FIG. 8: Diagrams representing 1)2.87^ 



Due to the factor \/Ecji}^) ~ l/i?ci(k) in the integrand, {Ncpix^)'''^^^'') A(p,x°) its magnitude much smaUer 



{c]cl) 



than that of (-/Vcp(a^°)^^')A(p,a:5)- 



From the model calculation explained above, though more detailed numerical analyses are necessary, we are possible 
to regard the magnitude of {Ncp{x^)'''^^) a(p,x°) dominant over other {Ncp{x^)''"^) Aip.x'l)^ 'i i-^- 



{Ncp{x°))a(p,x») - {]^Cpix°)''^^)A(p,x°)- 

If we add by hand the term due to the lifetime of A-particle in the way as 



(2.38) 



^zEa{p)z° _^ ^{^EA{p)-^^)z'> ^ 



(2.39) 



in H2.35|l . we obtain 

{Ncp{x'f'')Aip,x») = v\fABc.\'Y.KrK('^y'^' 



Jj 



(Pk 

(2^ 



^^{Ec,{\^)-Eclmt 



Ec,(y) Eci{^) 



2EB{q} 



gi(£;A(p)--EB(q)--Ecj(k)+irA(p)/2)t _ I ^i(-EA(p)+EB{ci)+EciCk.)+iTA(p)/2)t _ ^ 



Ea{v) - Eeiq) - Ecj{^) + *r^(p)/2 S^(p) - ^^(q) - Ecii^) - iTA{v)/2 



, (2.40) 



q=p-k 
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Under the condition TA{p)t ^ 1. we have 



{Nc,{x'f')Ai,^.o, - V\fABC.\' E KrKl'^l^Z^f J ^'^'^^^■^'^-'^^"'"'G^b.Kp, q = P - k, fc), (2.41) 



3,1 

1 1 



1 



2i?B(q) [Ea{p) - Eeici) - Ecjiy) + *r(p)/2] 



[Ea{p) - i?i3(q) - i?ci(k) - ir(p)/2] 



(2.42) 



E.a. p.40() should be compared with E.a. p.l5|) . It will be meaningful to examine whether Gab-jI depends weakly on 
and then, whether 12.41|l depends on t, characteristically different from the standard formula \i.b\ . 



III. EXPECTATION VALUES OF FLAVOR-NEUTRINO CURRENTS WITH RESPECT TO ONE-PION 

WAVE-PACKET STATE 

A. Starting Relations 

We examine the case where flavor neutrinos are produced through a boson decay such as 7r+ —t /i+ + i^^; for 
convenience, nuetrinos with only one flavor a are produced. As the initial 7r+-statc |5'((p) ;X, a;^)), we adopt, in 
accordance with Giunti et al.jlOj, 



|vl/((p) ;X,x?)) = jSpAliv, (p) ;X,a;?)-=l=4 (p) |0), 
J -^/^ii^^pj 

Alip, (p) ; X, X?) = A;(p, (p) )e-'P-x+'^^'(p)-? 



(V2^a03/2 



exp 



-^P^^-^p.X + ^i^.(p)x? 



(3.1) 



(3.2) 



Here, the operator a-„{'p) appears in the plane- wave expansion of the 7r^-field, as already given by H2.3|l . The state 
normarization is given by 



(vl/((p);X,:r?)|vI/((p)';X',:r?)) 

^2 



= exp 
x=x' 



f(X-X')^-— ((p) -(p)')'-* 



,,2 ,((p) +(p)')-(X-X') 



exp 



8<T 



■S;j((P>-(p)r 



(p>=(p>' 



1. 



(3.3) 



We have the following relations 
(i) 



(ii) 



Jd'pAlip, (p) ;X,x?)(-«^)A.(p,(p) ;X,x?) 
= Jd^pp\AAp,{p))\' = {p), 



d'pAlip, (p) ; X, x")(-z^)A,(p, (p) ; X, x°) 



(3.4) 



d^pE^{p)\AAp,(p))\' 
Jd^EAiP) ) + (v.) (P - (P) )]|A(p, (p) )p = K((p) 



(3.5) 



where 



dp 



p=(p) ^-((P>)' 
hereafter Ey^{{p) ) is written as {E.^) for simpHcity. 
(iii) By employing the above relations (|3.3|) - (|3.5|) and the energy and momentum operators of Tr^-field 



we easily obtain 



/ 2^ (^P^O ["JrP)"-(P) +/?l(-P)/5^(-P)]. 



(vI/((p);X,x?)| ff\ |*((p);X,a;?)) M^^^^ 



(iv) The wave packet in the coordinate space is defined by 



AliK, X?; X, xl (P> ) ^ / (^^^:(P' (P) ; X, x?)e^(P-^^(P)-") 



By using the approximate form of E-^{p) in 1)3.5(1 . we obtain 
A;(x,x°;X,xO,(p)) 
■ exp 



1 



where CTTrCTrx — 1/2. The normarization is 



z (p) • (x - X) - z (E^) (x° - x?) - fc-^^ " ^'^'^^ ^"^^ " "^"^1' 



4ct2 



which is the same as (|3.3p. 



(v) 



Jd'xAli^, x?; X, (p) )A,(x, x?; X', (p) ') 
= J dV;(p, (p) ; X, a;°)A,(p, (p) '; X', x?), 



y^^x X |A,(x, X?; X, X?, (p) )|2 ~ X + (v„) (x" - .T?). 



Next, the relevant effective weak interaction, causing 7r+ — > + i/g. decay is written as 

Hwix) = -Cint{x) = -[P„{x)J^{x) + J^{x)iy„{x)]; 
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(3.6) 



(3.7) 



(3.8) 



(3.9) 



(3.10) 



(3.11) 



(3.12) 



(3.13) 



the effective forms of Ja{x) and Ja{x) have been given by p.22|l and p.23|l . The spin 1/2-field V'(a^) is expanded in 
the interaction representation as 



hclicity r 



[a{kr)e'''''u{kr) + P'< {kr)e-'''''v{kr)] . 



(3.14) 



The concrete forms of u(kr) and v{kr) in Kramers representation^^ as well as other details are given in Refs.0 and 
. 

The expectation value of the flavor neutrino current j'^^^{x), given by (|1.26|) . when the initial state is 
|^^((p) ; X, xj)), given by H3.1|l . is expressed as, in accordance with (|1.19l) . 



(vl/^((p) ; X, x?) I 5-i(x°, X?) : jf^) (x) : ^(x°, x°) \^.i{p) ; X, x?)). 



(3.15) 



This consists of a sum of two parts in the lowest order of Ti.w ] 



f(a;°;(p),X,a;?)^i)^ 
f(a;°;(p),X,a;?)^„)^ 



X nwiy)nw{z) : : + : : 'Hw{z)nw{y)\ |vl',((p) ; x, 

In the following we examine H3.16|l . because this includes a dominant contribution. 



B. Form of E^j-^^ 
As shown in FIGEl £^Y)p consists of two parts expressed as 



■£(x,x";(p),X,a;?)- 
f(x,xO;(p),X,xO)« 



(i-i)p 

a 

(I-2)P. 



= (vI^,((p);X,a;0)| 



d'^zd y 



Ja{z)Va{z) ■■ 3(p){x) ■ V^{y)Ja{y) 

Va{z)Jo{z) ■■ i(p)(a;) : Ja{y)va{y)^ 



|vl/,((p);X,X?)) 




(1) ^(I-l)p 




We obtain 



FIG. 9: Diagrams representing 113.18^ 

£(x,a;";(p),X,a;°)^i_i)^ 

= / d^^A (*^((p);X,x?)|z/:,9,4(a;) 

JxO Jx° 

X {0\l{z)v''{0\ v„{z) : jl^{x) : ^^^(y) |0)z;'^/,(y) |0) • t^ddM^) |*.((p) ; X, x?))^ 



d^zd^y Id^p' Id^p 



=A.(p',(p) m;(p,(p)) 



(27r)V2-B.(p')2S.(p) 



^ g-»p'-(z-X)+i£;,(p')(z''-a:0)g»p-(y-X)-j£;,(p)(yO-an?) 
X ^'i'{'l+k'j)z-i{ki+q)y ^i{ki-kj)x 
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where ki = (k, iw/(k)), W((k) = ^/WTm^ ; 

= X! ^<^iis)^'il + j5)uj{k'r') {-i)uj{k'r')j''ui{kr) ■ ui{kr)^(l + ^5)vaiqs). 
By taking the new parameters 

and adding by hand the factor of tt hfetime in the same way as (12.3811 . we obtain 



£(xO;(p),X,:E?r 



(i-i)p 



dZ°dY° d^p' d^p 



(^3_^j3ygiz-(-p'+q+k')+iy-(p-q-k)giX-(p'-p)+ix(-k'+k) 
^ gjZO(iS„(p')-B<,(q)-c^j(k')+ir^(p')/2)giyO(-i5,(p)+_E„(q)+u.,(k)+»r,(p)/2) 



X e 



it(wj(k')— wi(k)) 



(3.20) 



(3.21) 



(3.22) 



The factors depending on and y", when we perform the Z^- and y '^-integrations and require the conditions 



ff^(p') > 1 , tV^{p) > 1, 



lead to the factor 



(3.23) 



£;,(p')--B.(q)-Cc>,(k')+i 



r.(p') 



r.(p) 



q=p'~k' 
=p k 



Then we see that possible candidates having a connection with the oscillation in R.H.S. of 13.22|) are 

i) e»("^(k')-c^,(k))t g»(p'-p).x+«(-k'+k).x ^ g«(-k'+k).L ^j^j^ L = X - X. 

We will examine the situation in the next section. 
Similar to f ("j.^-jp, we obtain 



£(xO;(p),X,a:?)^i_2)p 



Jo 



dZ°dY° d^p d^p' 



, A(p, (p) )a;(p', (p) ) 

(2^)V2i?.(p)2i?.(p') ' ' .VP,\P// 
$:z:fZpf*z;/^Z^/^*sg;(,;p,p^fc,fcO-e^<''-'''^-<^-^^e^^--('')+-^(^^ 



^ gizoc-E^CpO-B. (q)-(^j (k')+ir^ (p')/2) (p)+E„(q)+c^, (k)+ir, (p)/2) 



q— — p' — k— — p— k' 



(3.24) 



(3.25) 



(3.26) 



here, 



Bi';);(g;p,p';fc,fc') 

= X! '"<t((?s)/^(1 + 75)wi(fc?') ■ ivi{kr)'y°'Vj{k'r') ■ Vj{k'r') ^' {I + '^5)ua{qs). 



(3.27) 
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The factor corresponding to (|3.24(l is 



r.(p') 



E^{p) + E,{q) + uji{k) + 



r.(p) 



q=-p'-k' I 
= -p-k 



(3.28) 



Though detailed numerical analyses are nescessary, we can expect the main contribution in £(i^p comes from S^j -^^-^^. 
In the next subsection, we perform further evaluations of f ('j_x)p- 

C. Evaluation of S^i.^^p 

With the aim of examining the effect due to the wave packet of 7r+, we perform the integrations in R.H.S. of (|3.22() 
with respect to p' and p firstly; then the z and y integrations are carried out. 



1. p' - and p-integrations 

We consider the p integration appearing in (|3.22() : 

Jd'pAlip, (p) yy°i-EAP)+^rAP)m+^iy~^)■P ^(p)^ 

and take the approximation, similar to H3.5(l . 



g^F°(-£;,(p)+^^,(p)/2) _ ^^yO (- (BJ _ (v, ) ■ (p- (p))+i(r„ ) /2) . 



(3.29) 



(3.30) 



here, (T.^) = r„((p) ). Also F{p) in the integrand is assumed to be replaced approximately by F{{p) ); then we have 



- ^((p> ) 



= Fi{p)) 



27rcr^]3/2 



exp 



(P-(P))^ 
4a2 



*(p-(p))-(y-x-r" (v.)) 



+z(p) .(y-X)+*y"(-(K) +^'-^) 



27ra^]3/2 



exp 



4a? 



(p-2za2(y-X-r" {v^))f 



X exp 

~ F((p))(2a^\/2^)3/2 exp 
Thus, H3.22|l is written as 



a^(y - X - y° (v.) f + I (p> • (y - X) + ^Y'{- {E^) + ^^^) 



-al{y - X - r" (v.) f +^{p) ■{y-X)+ rY\- {E^) + ^^) 



£:(x,x°;(p),X,x?)^i_i)^ 

i y ^2^)32(£;.) 

,1/2^1/2* 1/2 „l/2*(l)a 



d-'q d-'k' d*k 



E ^^f ^rXr^r i^Sri-z; (p) , (p) ; fc', ma.v2^\ 



X giz-(-{p)+q+k')+iy((p)-q-k) 

X g«z°((-Bx)-B.(q)-c^,(k')+»{r,)/2)g*y°(-{£;,)+iJ„(q)+c^,(k)+j(r,)/2) 



X e 



ix-(-k'+k)+it(Ljj(k')-w;(k)) 



exp [-a^(z - X - (v.) )2 - ^^(y - X - y° (v^) )2] . 



We see 



y - X - r" (v,) =0, z-X-Z° (v„> = 
correspond to the classical relations as shown in FigllUI 



(3.31) 



(3.32) 
(3.33) 
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(X,t = 0) 



FIG. 10: Classical Relation iT^ . 



2. z- and y -integrations 



By employing 



-al{y - X-rO (v.) )2 + zy . ((p) - q - k) 

= -aliiy - X - r" K) ) - ^((p) - q - k)]^ 



4aE 



((p> -q-k))2 + i(rO(vO +X)((p) -q-k), 



we obtain from H3.32|l 



J(x,i°i(p),X,i»)J.,,, 



xi?W-(g; (p),(p>;fc^fc)[2a,^/2^] 



X exp 



-^((p)-q-k)f-^((p)-q-k')f 



(3.34) 



^ g»Z°({i5,>-B„(q)-^,(kO+i(r,)/2+{v,)-(-{p>+q+kO)giy'\-(S,> + S„(q)+^,(k)+i(r,)/2+(v,> 
y ^l(x-X)-(-k'+k)+it(i^i(k')-"!(k)) 



(3.35) 



5. k- and \l' -integrations 

Similarly we can perform k- and k'-integrations approximately. For a fixed K, which is defined by 

K = (p> - q 

the main contributions to the k and k'-integrations come from the regions 

k ~ K, k' ~ K. 

Then, we are allowed to employ the relation 

duji 



9k 



(k-K) =c^z(K)+vz(K). (k-K), 



(3.36) 
(3.37) 

(3.38) 



K 
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and that for ^^(k'). Then we rewrite the exponent of H3.35|l as 



- -^(k - Kf - iioi(k){t - r°) - iy" (v,> • (k - K) + ik • (x - X) 
-^(k - Kf + t{k - K) . [L - v,(t - yo) - yo (vO ] 

- iu;i{K){t - Y°) + iK ■ L 
= [(k - K) - 2talFi{K; L, t, Y")] ' - a^FfiK: L, t, V') 

-iwi(K)(t-y°)+iK-L, 

where L = x — X and 

F;(K; L, t, y") = L - v,(K)(i - - (v.) . 
By employing l|3.39|l . we obtain from 13.35|l 



£:(x,a:0;(p),X,x?)^i_i)^ 



JO 



xBl^l'^iq; (p),(p);K,K)[2a,V2^]' 



(271)122 (ii;^ 



1/2* ^1/2 ^1/2 7I/2* 

pj '^pi ^<yj ^<Ti 



^ ^lZ"((E„)-E,(c^)+^{T,)/2)^lY•■>(-{E,)+E,(c^)+^{T,)/2) 

X exp [-al¥, (K; L, t, - a^F;(K; L, y°)2] 

iwj(K)(t-Z'')-iw;(K)(t-y°) 



X e 



K={p>-q 



(3.39) 



(3.40) 



(3.41) 



In R.H.S. of (|3.41() . the L-dependent osciUation factor in H3.35|l . exp[i(k — k') • L], disappears. The L-dependence 
is expected to be recovered due to F^- and F^-terms through the time integration in (|3.41|l . In order to see this 
possibiUty quaUtatively, we evaluate roughly the Y^- and Z^- integration in (|3.41(l under the classical relations (as 
shown in FIG El) 



F,(K;L,t,yO) = F,(K;L,t,Z°) =0 



(3.42) 



L('=X-X) 



(x,t) (y) (X,t = 0) 



FIG. 11: Classical Relations 



Under the parallel condition, where the relevant momenta, (p) and (K) = (p) — q, are parallel to L, the solutions 
Y^ and Z° of (|3.42|l are given by 



yO^jntL . zl^Y^{v,^v,), (3.43) 



t-Yy = 



t-Z^,=t-Y^{vi^v,). (3.44) 
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By assuming < {Y^, Z^} < t, and by taking the rough evaluation of integrations 

f dZ° fdY'^---^ fdZ° fdY%{Z°-Z^^)S{Y°-Y,°), 
Jo Jo Jo Jo 



we obtain 



arm 



i2a.f E Jd'lBl^itiT, (P) , (P) ; K, K) 



exp [(7'jj)clas] 



K=(p>-q 



where 



(T'.Ocias ^ zK(K) - c^mt - Z"A^SE, + ^) + Y,°{^SE, ^-^), 



with 5Ej = i?o-(q) + ^j(K) ~ (i?^) . We obtain easily 



Im{Vji)c\as = (i - (v,r) t) 



V, - VI 



[Vj - (v^) ){vi - (v^) ) 



VjVl 



K + (E^) - E„ 



(3.45) 



(3.46) 



(3.47) 



This unfamiliar form suggests a more precise evaluation of the integrations in (|3.41() is necessary, which we will 
consider as follows 



4. Z"- and Y^^ -integrations 

For simphcity, we write Fj(K; L, t, Z°) as Fj(Z°), and define Z^ by 

d . „ , 9I 

= 0. 



J.(F,(Z°))" 



Then we obtain 



(vj - (v^) )■ 



:(v,-(v.))-(v,i-L) 



= 2a2^.a2(v, -(v,>).(v,t-L), 
a^[F,(Z")]^ = - Zy + a^[F,(Z°)]^ 



where 



(3.48) 



(3.49) 
(3.50) 



a^(v,-(v.))^ = 



2<' 



()3.41|l is rewritten as 



£{L, t- (p) )^i_i)^ ^ £(x, x", (p) , X, x?)^i_i)^ 



M^Jd'^T. ^i^r^Mr^;^, (P> , (P> ; K, K) 



X exp 



rO\1 



i{ujj (K) ~ a>,(K))t - d^F, (K; L, t, Z")' - (7^F,(K; L, t, Y'^) 



7i ,r/.. i^On 



-^{iSEj 



-Y\^5E,^^^) 



X M,(i;Z") •M,(i;r")*, 



(3.51) 



(3.52) 
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where 



Mj{t;Z") = / dZ"cxp 



2<- 



The integration form of (|3.5^^(l is expressed as 



M,(i,^'0 = V2<exp 



where 



(3.53) 



(3.54) 



Jj{uj{t);aj, l3j) = / d-Uj exp[— (m^ + ia + /3j)^], 
Jo 

with Uj = Z"/(\/2cr7rj) and Uj{t) = t/{\/2aTrj)- We define the two functions as 



(3.55) 



fc{y;a,P)^ / (iuexp[a^ - (li + /3)^] cos(2q!(u + /?)), 
fs{y;a,f3)= / duexp[a2 - (u + /3)2] sin(2Q(u + /?)); 



then, (|3.54|) is expressed as 



where 

0(?/; a, (3) = arctan 
Thus, H3.52|l is expressed as 



I -^i ("j (0 ; > /^j ) I exp [-iQ{uj{t)- aj , /3j )] , 



/c(y;a,/3) 



(3.56) 



(3.57) 



(3.58) 



X Bl]l'^{q; (p) , (p) ; K, K) (2a,, a,,) exp 
X \ Jj{uj{t);aj,Pj)\\Ji{ui{t);ai,f3i)\ exp 



<.(z<5i?, 



(3.59) 



K={p)-q 



where 



[F,(^°)]' = 



Z"i^SE,+ ^^)-n-^SE, , ^ „ 
^ [(v, - (v.) nv,t - L)2 - {(v, ~ (v.) ) . (v,t - L)}2] 



(vj - (v.) )^ 



(3.60) 
(3.61) 



with = K/^K2 + m2, K = (p) - q. 

Though the t- and L-dependent part H3.59|l is comphcated, it becomes a httle simpler when the parallel condition 
is assumed. From [Fj{Z_^)^ + Fi(y")2]para = 0, which leads to the classical relation 



[Z\ 



jvjt - L) 
Vj -M' 



(3.62) 
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we obtain Im['Pji]para, 

Im['Pji]para = {L - (u^r) t) 



Vi + Vi - iV^ 



VjVi 



(3.63) 



which is the same as (|3.47|) . 

We can examine the quaUtative structure of H3.59|l under the condition of neutrino high energies in comparison 
with m^ 's. This is to be explained in the last section Sec.V. 

IV. EXPECTATION VALUE OF NEUTRINO CURRENT WITH AN ADDITIONAL OBSERVED 

PARTICLE 

As pointed out in Sec. Ill, it is worthwhile for us to examine the structure of the expectation value of the flavor 
current, where the charged lepton observation is taken into consideration. 



A. Boson model; the plane wave case 

We first consider the boson model in the plane wave case with the aim of seeing the way how to take into account 
the observation of an additional particle. We examine the same model as that in Sec.IIB ; 



Wint = -fABCa(l)%„{x)(j)Ba{x)(f)A{x) + B.C. 

'Hii^) = ^lJx)Mx), H2ix) - n\ix) - 4{x)^cAx), 

From the expectation value H2.28|l Cwe pick up the partC corresponding to FIGEfl), i.e. 

£{A{x'}y,Cp-cm{x)t'^ ^ d^z A (0| aA(p)i/2(^) : j'f'Hx) : H,{x)al{p) |0)^ 



(4.1) 
(4.2) 
(4.3) 



d^z / A (0| aA(p)4(^) • 1 • J.(y)a!i(p) |0) 



x(O|0c.(^):jf''^(x):0L(?/)|O)^ 



(4.4) 



The first part in the last integrand includes the propagator (0| <l>B„(z)4>Ba{y) |0). This can be effectively set equal to 



U'LW'^B.(y)|O) = (O|0^j,(z) 



d^q 



2EBici) 



/3L(q)|O>(O|/5B.(q)0B.(2/)|O) 



nBa{q)d^q 



<i>Ba{y)W 



(4.5) 



with -UBaiq) = l3^B^{q)(3Ba{q) / {2EB{q))Cas seen from 



\PBaiQ') 



-jnBa{q)d\ /3tj^(Q) |0)™" = J-^-2Eb{Q')S{Q' - q)2EB{Q)S{Q- q) 



I/3b.(Q')/3L(Q)|0)- 



£{A{x°y,Cp-Cm{x),BM)-ohs)i'~^^ = / / d^zd%{0\aA{p)4{z){-nBAq))JAyWA{p)\0r 

JxO Jx° 

x(O|0c.(^):jP^(^):0L(?/)|O)™^ 



^B.vv«;,w (4.6) 
ThusCwe can regard the expectation value of the flavor current with the additional particle observation is given by 



(4.7) 
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(See FIGd) 

The concrete form of (|4.7|l is given by 



(27r)62wj(k')2w,(k) 

* ''*d^0^yO_^4BC^ ^^1/2^1/2*^1/2^1/2* ^.(., (k) -.Hk))t 



{2^f2EB{ci) 



2wj(k)2wi(k) 



2k 

z(a;j(k)+c^,(k)) 



^l{EA{p)-EB(.q.)-uj,{k))z°-iiEA{p)-EB{q.)-i^,{k))y° 



k=p q 



where t = x'^ — x'j. We see ()4.8|) leads certainly to 

fd^x jd^k i£{A{x°j);Cp-Cm{x), B^i^yohs) 



4 ~ 



^Cp(x'')W)a(p),.o 



(4.8) 



(4.9) 



given by 

Due to the plane- wave calculation, the x-dependence disappears in 14.8|l : thus, there is no spacial coordinates 
specifying the observation. In the next subsection we examine the model calculation by employing the wave packets 
for the external particles ( i.e. for 7r+ and the "observed" i^CT-particle). 



A 




riBaiq) 




FIG. 12: Diagram reprenting 14.71 



A 



B. Boson model; the wave-packet case 



First we have to take out the part, which is suitable for the present purpose, of 



(*^((p);X^,a:?)| 



d^zd%Hint{z) : jf'Hx) : Hint |*a((p) ; X^, x?)) 



\ \ con 



Here, we choose the wave-packet state in the same way as H3.1|l : 

|vl/^((p) ;X^,x?)) - /d3pA^(p,(p) ;X^,x?)-=l==aIi(P) |0), 



A^(p,(p);X^,x?) = 



■ exp 



^5— M]!_,p.x^ + *i?^(p)x? 



(4.10) 

(4.11) 
(4.12) 



(V2^(7a)3/2 

This state is normalized in the way as noted by (|3.3I) . (|4.11|) is regarded as the state obtained through substituting the 
wave-packet state for the plane- wave state . ctti(p) |0). In the same way, the wave-packet of Be incorporated 

in H4.5|l through the substitution 



2i?B((q)) 



/3L((q> ) |0)(0|/3s.((q) ) |*B((q) ;Xs,x"s))(*B((q) 



(4.13) 



23 



with 



|vI/B((q);XB,4)) = jd\A*s{ci, (q);XB,x°) 



|0), 



(4.14) 



where A*^ is defined similarly to H4.12|l . 
Thus, we are led to examine 



f [A((p) ; X^, x?); Cp-Cur(x), B,((q) ; ^b,x%)]. 



Nia-i) 

d4zd4y|/ABCaP(*A((p) ; X^, |vl/^((p) ; X^, a:?))™'^ 

X (0| W I*i3((q> ; Xs, xO,))(vl/s((q) ; Xs, xO,)| cj^B^v) |0) 
x(O|0c.W:jP^(:^):0L(y)|O)™^ 



(4.15) 



We add a remark on the choice of ()4.13f) . As noted above the state normalization is changed from (0| a^(p')at(p) |0) 
2Sa(p)<5(p'-p) to 



(vl/^((p') ;X^,a;0)|*A((p) ;X^,a;0)) - exp 
In the same way, we adopt (14.13(1 which satisfies 



(4.16) 



(0|/3b.(Q') d?{q) |vl/B,((q);XB,x^))(*B.((q);XB,4)|/3L(Q)|0) 



dP {q) V2i?i3(Q')2i?B(Q)A|j(Q', (q) )Ab(Q, (q) ) 



v/2£;s(Q')2-Bi3(Q)exp 



8-1 



instead of adopting 



(I4.18|) satisfies 



AAJj(q,(q) )AB(q,(q)) /^^(q) |0) (0| /3s.(q). 



(0| PBa{Q') (^j (q) 6IHj) /3L(Q) |0) = 2£;b(Q)5(Q' - Q), 

which corresponds to the plane-wave normalization. 
We obtain from 1(4. 15() 

£[A{{p) ; Xa, xf)-Cp-GuT{x), 5,((q) ; Xb, 4)]['-^) 
= J^^ J^^ d^zd^y\fABCa\^ /"^'^'/"^'^ jd\' jd\ jSk' jd^k — 



(4.17) 



(4.18) 



(4.19) 



(27r)V2£^A(p')2^A(p)2^s(q')2-BB(q) 

5^ zl'^'zlfzfzlfAAip', (p) )Al(p, (p) )AJ,(q', (q) )AB(q, (q) ) 



gi(-p'+q'+k')-z+i(p-q-k)-ygi(p'-p)-XA-i(q'-q)-XB+i(-k'+k)-x 



(27r)62wj(k')2w,(k) 

X exp [i{Ea{y>') - Eb{c{) - ^,(k'))(z" - a;?) + «(-£;a(p) + ^^^(q) + ^Kk))(j/" - a;?) 
+*(c.,(k') - u;,;(k))(a;° - .t?) + *(£;B(q') - £;B(q))(4 - 2:?)] • 



(4.20) 
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We perform the integrations with respect to p', p, q' and q by empfoying the approximate relations as already used 
in Sec. Ill fe.g. l|3.29|l ). and obtain 

£[A{{p) ; Xa, x% Cp-Cm{x), B,((q) ; Xs, x^)]!'"^) 



dZ'dY" I Sz I cPy I Sk' / cPk I/. - - 



oJo J J J J {2TTf2{EA)2{EB) 



ABC<j\ 



k'+k)-x+j:(cjj(k')-Wi(k))t 



(27r)62wj(k')2w,(k) 

X exp [i (p) • (y - X^) + {-I (Ea) - {Ta) /2)Y" + t {p) ■ (-z + X^) + (^ {Ea) - {Ta) /2)Z° 
+ ■ i-y + Xb) + * {Eb) (r" - ) + z (q) • (z - Xs) - * (iJ^) (Z" - X^) 
-/^(y, r") - h{z, Z") + t{^Lu,{k')Z° + Lu,{k)Y" + • z - k • y))] , (4.21) 

where the weak-decay width of A-particle is added by hand; t — — x'j — ~ x'j , Z^ = z'^ — x^ , Xg — x^ ~ x^; 

h{y, y") = <j\{y -Xa- {va) Y^ + a|(y - X^ - (v^) {Y^ - X%))\ 

With the aim of execute the z- and y- integrations, we first examine the structure of 

exp [-/i(z, Z") + i{- (p> + (q> + k') • z] , 

which is included in the integrand of 14.2111 . The concrete form of z , defined by 

5/i(z,Z0) 



9z 



0, 



is given by 



then, we obtain 



[a\{XA + Wa) Z") + oHXb (vb) {X% Z'))] 



(4.22) 
(4.23) 



(4.24) 



(4.25) 



Mz,Z") = 



'A-B [Y^B^z\X%)f 



(4.26) 



with 



and 



F^^(ZO, ) =Xb- (vb> {Xl - - (VA> ^" - X^, 



Mz,Z") = (a2 +4)(z-z)2 + /j(z,Z°). 
We see exp[— /i(z, Z'^)] is maximum at z = z, and 

/i(z, Z°) = ^ ¥^^{Z°, X%) ^ 0. 



F"^-^(Z",X^) leads to the relation among classical trajectories. 



By noting 



exp [-{(t\ + a|)(z - z)2 + i(~ (p) + (q) + k) • z] 
exp y^ji^ I ^ -y, ^ 2(^-2 ^ ^2 
(-(p) +(q) +k)2 



-K + 4)<^z-z- 2^ (-(p) +(q> +k) 



«(-(p) +(q) +k) -z 



(4.27) 



(4.28) 



(4.29) 



(4.30) 
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we perform the z and y integrations of (|4.21() and obtain 

£[A{{p) ; X^, x?); Cp-Cur(x), S,((q) ; Xb, 4)]!"''^ 

OJO J J {2nf2{EA)2{EB) 

E ^1/2^1/2* 1/2^1/2* Z^(fcj- + ki)b i(-k'+k)-x+i(a;j(k')-i^,(k))t 

« '^^ (27r)62u;j(k')2u;,(k) 



TT 



X exp 



- /j(z, Z°) + (p) + (q) + k') • z - h{y, + z((p) - (q) - k) • y 
1 



{(-(P> +(q> +k'f + (-(p> +(q) +k)2} 



1 



+ -(i^s) -c.,(k')) + *r°(-(i?^) +{Eb) +c.,(k))--(r^)(zo + y") 

Next we proceed to integrate with respect to k' and k. By employing the definitions 

(k) ^ (p) - (q) , {lo,) ^ ), (v,) ^ v,((k) ), 



and 



exp 



h{x, ZO) + i{- (k) + k') • z + iZ'^i iEA) - (Eb) - {uj) J - (v,) • (k' - (k) )) 

1 



exp 



ik' . X + ^{ + (v,) . (k' - (k) ))t - \ , (- (k> + k')2 



- h{z, Z°) - 



1 



{k' - (k) + 2t{a\ + 4) . (-Z + (v,) + X - (v,) t)y 



(ai + al) {-z + (v,) + X - (v,) t}' + ^Z\{EA) - {Eb) ~ (^,) ) 
J (k) • X + i {ujj) t , 



we obtain from H4.31|l 



£[A{{p) ■ Xa, X?); C,-Cur(a;), B,((q) ; XB,x%)]i'''^ 



7I/27I/2* l/2„l/2* 



"'0 



(2^)62 (Ea) 2 {Eb) 



(Ja + '^b 



x[2Jn{al + al)]' xe 



2 ^l6 ^ i({i^j>-(wi))t 



X exp 



Gj{Z°, Xa; X^.Xb; t, x) - Gi(F°, X^; Xg,XB;t, x) 



+ *Z-'((i?^> - {Eb) - {uj) ) + iY^{- {Ea) + {Eb) + {to,) ) - - (r^) (Z" + Y") 



2(k) 

+ {uj,) ) 



(27r)62 (wj) 2 (w,) 



where 



(ZO, X^; X^, Xs; x) = h{z, Z") + {a\ + 4){x - z - (v,) {t - 

In order to see the properties of the integrand, it seems convenient for us to define 

Ff{Z') = X - (v^> ZO - (v,> {t - Z") - Xa. 
Ff^(Z",X^) =F^^(Z",XO) -F/C'(^") 

^Xb- (vb) {Xl - Z") + (v,) {t - Z") - X, 



(4.31) 



(4.32) 



(4.33) 



(4.34) 

(4.35) 

(4.36) 
(4.37) 
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where F^^{Z'^, X'^) is given by (|4.27() . The factor x - z - (v^) (t - Z°) inchided in Gj, given by (|4.35(l . is expressed 



as 



X - z - (v,) (< - ZO) = F/^(Z") + (v^) Z" + 
1 



B L 



due to 1(1:^ 



thus we obtain from l|4.26|l and H4.35|l 



(4.38) 



(4.39) 



Here we give a remark on the possibihty of the minimum of Gj to be 0. The condition Gj = for Z'-^ = Z^^^^ is 



Fr(^cias) = 0, and Ff^(ZlJ = 0, [due to ^] 



leading respectively to 



1 



clas 



clas 



:((v,> -(v^>).(-x + XA + (vj)t), 



Jf^-y^J^ii^i) ~ i^A) ) • (-X + Xs - (vs) X% + (v,) t), 



(4.40) 

(4.41) 
(4.42) 



together with the componets of — x + X^i + (vj) t and — x + X^ — (v^) Xg + (v^) t perpendicular respectively to 
(vj) - (va) and (vj) - (vs) to be equal 0. Eas. ljOUjl means F^^(2'°i^s, ^s) = [due to (|07jl ]. which leads to 



■^clas 



((va) -(vb))^ 



,((va) - (vb) ) • (Xb - (vs) Xl - Xa), 



(4.43) 



together with the component of X^ — (v^) X^ — Xyi perpendicular to (va) — (v^) to be 0. These Z^^^^^s give three 
kinds of the classical relations. We see, however, these Z^^^^^s should not be taken to be equal especially due to the 
j-independence of 14.43|l . 

It is necessary for us to consider how to perform the Y'^ - and Z*^-integrations over a finite time-interval; this is to 
be considered in the next subsectionCwhere the minimum of Gj is examined. 



C. Time integrations 

We examine in detail the structures of the integrand in (|4.34f) . 



1. The minimum point Z^^ of Gj 
We write for simplicity Gj as Gj{Z^) and rewite it in the form 



G,(Z°) = a,(Z°-ZY + G,(^°). 



(4.44) 



Z^ is determined by 



= = 2 



al((v,) - (va) ) • Ff ^ {Z^) - 4((v,) - (vb) ) • (v^.) f ^ (Z», X%) 



due to 



{4((v,> -(v^>)^ + 4((v,> -(vB>f}^° 



which leads to 



4((v,> - (v^) ) • Ff^(O) + 4((v,) - (vb) ) • Ff^(0,X^) 



here 



with 



11 1 



By noting Uj to be given by \{(PGj/d{Z^Y) = l/(2cr|), we have 



2a| 



^. Another expression of Gj{Z^) 
We write the points Z^'s giving the minima of (F^*-^)^ and (F^*^)^ as T^^^ and , respectively 



We obtain 



then H4.46|l leads to 



By noting 



FfiTf^) ■ ((v,) - (v^) ) ^ 0, Ff^{Tf^, Xl) ■ ((v,) - (vb) ) = 0. 



AC 



1 



rpBC _ 



1 



((v,) -(v^))-(-x + X^ + (v,)t) 
M'^m((v,> -(v^>)-(-F/^(^° = 0)), 

((v,) - (vs) ) • (-X + (v,) t - (vs) + Xs) 



T 



AC 



T 



BC 



^2 rpAC , ^2 rrnBC 



1 9^ 

2 9(Z0)2 







r 1 -| 






(pBC)2_ 
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we have from (|4.39() 



It is easy for us to confirm the equality 



1 



1 



2af -' 2ia^^+a^sj 



2a 



because, due to (|4.52|l we have 

L.H.s. of g:^ 



1 



Thus we obtain 



1 



+ 

R.H.S. of 



n 2 



rpSC\ 



(4.54) 



(4.55) 



(4.56) 



(4.57) 



3. The expression of 



H4.34|l is expressed as 



£[A{{v) ; X^, X?); C,-Cur(a;), 5,((q) ; X^, a;0,)]['-^' 



2 (i;^) 2 (i;^) 2 (cc;,) 2 (cc-,) 



2(k) 



X exp 



- Z°iiAE, + ^)+Y''{iAE, + ^) 
X Ij{x,t;XA,XB, Xg) I*{y:,t;XA,XB, Xg), 



/j(x, t; Xa,Xb, Xg, 



X" 



dZ° exp 



(Z0-Z ")2 



(z° - z°)(iA£;j 



(Fa) 



here, we 

AE.J = (Eg) + {lo,) - (Ea) . 
and required a physical condition on the time ordering 

x"^ >x%>{y°,z°), i.e. t>X°g>(Y°,Z'^) 



In the same way as Mj{t;Z^) given by (|3.54|l . we can write Ij as 



Ij = \/ 2(7^ exp 



1 2^- A TP , (^a) ^2 



J,(7/,(X°);fi„/3,), 



(4.58) 
(4.59) 

(4.60) 
(4.61) 

(4.62) 
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where a, = (TjAEj/V2, j3j = (-Z" + icr^ )/{^/2cfj) and Uj{X%) = t/{^/2(j^j). Then, is expressed as 



£[A{{^) ■ y.A,x'\)-Cp-Gnv{x),B„{{q) ; Xs, x^)], 



m(/-i) 



E7l/2„l/2*7l/2 7l/2* I . |2/„2„2n 



3/2 



2 2 2 (wj) 2 (w,) 



2(k) 



exp 



X \ jj{uj{X%);aj,(3j)\ | Ji(ui(X^); 5^, /3i)| exp 



*(e(^.,(X°);5„/3,)-e(j^*)) 



(4.63) 



where 

Q,, = i{{uj,) - {cj,) )t - Gj{7P) - G,(rO) - Zt\%AE, + il^) + rO(zA£;, - ^). (4.64) 

We see (14.64(1 has structures similar to those of (|3.59|) . and examine in Sec.V its quaUtatives oft- and L-dependence 
under the condition of neutrino high energies. Here we add a remark on Gj{Z^) as follows. From H4.48|l and l|4.57|l . 
GjiZ'^) is given by 



(rpAC rpBC\2 



(4.65) 



As noted in the subsection B of Sec. IV, we can not take two (Z°)cias's (given by (|4.41|) and (|4.42|) 'l to be equal. This 
means Tj^'~^ and T^'^' , given by (|4.5U|) and 14.51II . are not equal to each other, since the expressions of (|4.5(J|I and 
(|4.51|) are equal respectively to (|4.41|) and H4.42|l . Thus, due to (|4.52|l . we have 



3 ^ — 

Under the parallel condition we obtain 



3 ^ — 



{Tj )para 



y-^3 



{^3)t-L 



L ~ X — Xa , 



-{{^j)t-L + LBA-{vB)Xl), 



Lba — Xb — Xa, 



(ZO 



/para o , ^2 l"3B\^] 



IF^jB (^Z*" )para + (^^j A (Tj^*" )para] ; 



(4.66) 

(4.67) 
(4.68) 
(4.69) 



the minimum of Gj{Z^) is given by 



G'j((Z°)para) — 



;para V-'j ^paraj 



(4.70) 



DISCUSSIONS AND FINAL REMARKS 



We have investigated the structures of the expectation values of the flavor current with respect to an one-particle 
state which plays a role of the neutrino source. The present approach is certainly simple one to the neutrino oscillation 
without recourse to any one-particle state of flavor neutrino. In Sec. Ill, the expectation value of the flavor neutrino 
current with respect to one tt^ wave-packet state was given as (|3.59|l . In Sec. IV, we have examined in the boson 
model the same type of the expectation value with an additional "observed" charged particle, and shown that the 
main part of the expectation value of (|4.63l) has the structures similar to those of (|3.59|) derived in Sec. III. 

The expectation values H3.59|l and l|4.63|l have somewhat complicated structure in comparison with the standard 
oscillation formula; therefore, in order to make clear characteristic features included in (|3.59() and H4.63|) . we are 
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necessary to perform numerical calculations, which will be done in a subsequent paper p^ . Here we give some 
remarks in the following. 

It is worthwhile for us to examine the features of contributions from high-energy neutrinos to H3.59|l and 14.63(1 . 
First we considered the structures of the integration of the integrand of (|3.59|) in the case of m^/fc^ ^ 1 for all j's 

and 1 - (v^) ~ 0(1). 

As to the phase parts, we obtain from 



Im(7',,) para 

with {5E) = E„ + K - (E^) , and also 

-Q{uj{t);aj,/3j) + e{u,{t); ai, Pi) 
By using th j-independent quantity 



Ami L~t{v^) 



2k 1 



m 



k{l - (v^) ) 



(5.1) 



dVj_ 

dm'j 



dQ{uj(t)]aj, I3j) 



m^j + ij ^ *)• 



(5.2) 







' dQ{uj{t);aj,j3j)' 




dv 







mj—0 



\Mu{t);a,P)\ 



^/{a''f + (/3')2 sin(2a/3 + 9 - (^^) 



(5.3) 



where primes (') denote the derivative with respect to neutrino velocity, v, and 



we obtain 



1 I ^ \ 2 

Lo — arctan — , <P = arctan , , , , 



99 



dv 



(5.4) 



This is the additional phase coming from the finite-time integrals in H3.52|l . One can see from ((5.3(1 that this additional 
phase will vanish when (3 and u{t)+(3 become much larger than 1 as well as a. This result can be understood because, 
since /3's represent the ratio of the distances {Z^ ^t — Z^) between the peak of gaussian (Z°) and both limits of integral 
region (0,t) to the gaussian width (cTrj), /3's ^ 1 means that the gaussian is so sharp that one can regard the both 
limits as infinite. Then the finite-time integral is reduced to the real function. However, the situation will be different 
when /3's are smaller than 1 and a. Numerical calculations are expected to make a behavior of this phase clearer. 
As mentioned before, the condition 



k^ 1 



«1, 



(5.5) 



is necessary for deriving ((5.1(1 . By taking (v-^) ^ and {6E) /k ~ 0, we obtain 



2k 



(5.6) 



which is the same as the standard one, apart from K — K(q) = (p) — q. It will be instructive to evaluate typical 
(w^) -values a shown in TableUfa); by the way, rough values of K by assuming \SEj\/ (E.^) < 0(1/10) are given in 
Table Cfb). 



It is possible for us to extract information on the coherence length from \Jj{uj{t);aj, f3j) \ ■ \Ji{ui{t); a^, /3,;)| as follows 
as KeVji- This is one of the remaining tasks. 
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(E) (MeV) 


200 


300 


500 


1000 


(a) 


(p) (MeV) 


142.8 


265.3 


480.0 


990.1 


M /c 


0.71 


0.88 


0.96 


0.99 


(b) 


k (MeV)* 


74 


122 


211 


427 



* Values of ^{ml - ml)/{{E^) - (p) ) 
TABLE L (a) and (b) 



Next we give a remark on the structures of H4.64|l . In this case, (k) (= (p) — (q) ) is fixed when (p) and (q) are 
given. We have the phase parts in H4.63|l expressed as 

Im(Q,,)para = [((t^,) ~ (^.) )t - ^E, +Y^^E,]^^^^, (5.7) 
[e(u,(t);5j,/3j)-e(j^z)]para. (5.8) 

In comparison with the case (|4.t)4|) . some comphcation appears due to the difference between (^*')para given by (|4.t)9|) 
and (^")para given by H3.62|l . This difference arises from H4.b6|l . i.e. from T^'-^ ^ "^f^^ which leads to nonvanishing 
[Gj{Z_°) + Gi(y°)]para in R-e[Qj.i]pai.a and also to the additional -dependence of Z° through (t|^ and cr^^. 



[1] M. Buete, Oscillations of nutrinos and mesons in quantum field theory, [hep-th/OlOgiigl Sep, 2001; Phys. Rev. D66,013003 

(2002); Phys. Kept. 375, 105 (2003). 
[2] S. M. Bilenky and B. Pontecorvo, Phys. Lett. 61B, 248(1976); Phys. Kept. 41, 225(1978). See also the recent articles, 
[3] E.g., M. Fukugita and T. Yanagida, Physics and Astrophysics of Nuetrinos, edited by M. Fukugita and A. Suzuki (Springer- 

Verlag, Tokyo, 1994), p.l ; S. M. Bilenky, Lectures given at the 1999 European School of High Energy Physics, Cast 

Papiemicka, Slovakia, Aug. 22-Sep.4, 1999 

[4] E.g., C. Giunti, C. W. Kim, and U. W. Lee, Phys. Rev. D45,2414(1992); S. M. Bilenky and C. Giunti, |iiep^ph/0102320| 

2001 

[5] M. Blasone and G. Vitiello, Annals of Phys. 244, 283(1995); 245,363(E)(1995); Phys. Rev. D60, 111302(1999). 

[6] K. Fuju, C. Habe, and T. Yabuki, Phys. Rev. D59, 113003(1999); 60, 099903(E) (1999). 

[7] K. Fujii, C. Habe, and T. Yabuki, Phys. Rev. D64, 013011-1(2001). 

[8] K. Fujn, C. Habe, and M. Blasone, hep-ph/0212076 Dec. 2002. 

[9] C. Giunti, C. W. Kim, J. A. Lee, and U. W. Lee, Phys. Rev. D48, 4310(1993). 
[10] W. Grimus, P. Stockinger, and S. Mohanty, Phys. Rev. D59, 013011(1999). 
[11] T. Yabuki and K. Ishikawa, Prog. Theor. Phys. 108, 347(2002). 

[12] B. Ancochea, A. Bramon, R. Munoz-Tapia, and M. Nowakowski, Phys. Lett. B389, 149(1996). 
[13] M. Zralek, Acto Phys. Polon. B29, 3925(1998). 

[14] M. Blasone, P. Jizba, and G. Vitiello, Phys. Lett. B517,471(2001); M. Blasone, P. P. Pacheco, and H. W. C. Tseung, 
hep-ph/0212402, Dec.2002. 

[15] H. Umezawa, Soryushi-ron (Misuzu Shyobo, Tokyo 1953), ch.l3; Quantum Field Theory (North-Holland Pub. Co, Amster- 
dam, 1956), ch.l3. 
[16] E.g. E. Cornaldesi, Nucl. Phys. 7, 305(1958). 

[17] E.g. V. de. Alfaro, S. Fubini, G. Furlan, and C. Rossetti, Currents in Hadron Physics (North-Holland Pub. Co, Amster- 
dam. London; American Elsevier Pub. Co, Inc., New York, 1973), ch.2. 
[18] T. Shimomura, in preparation. 



